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Two flavour massless QCD has a second order chiral transition which has been argued to belong to the uni¬ 
versality class of the 3d 0(4) spin model. The arguments have been questioned recently, and the transition was 
claimed to be mean field behaved. We discuss this issue at the example of the 3d Gross-Neveu model. A solution 
is obtained by applying various well established analytical methods. 


1. CHIRAL QCD 

Massless QCD has a flavour singlet, chiral sym¬ 
metry that is broken spontaneously, but becomes 
restored at sufflcently high temperature. For 
two flavours the symmetry restoring transition is 
widely accepted to be of second order. According 
to arguments first given by Pisarski and Wilczek, 
the universality class belongs to that of the 0(4)- 
spin model in 3d. Beyond others, two main as¬ 
sumptions made are 

• It is the symmetry pattern that determines 
the universality class. 

• Dimensional reduction is complete at the 
phase transition, that is, quantum fluctu¬ 
ations along the temperature torus do not 
renormalize the classical 3d IR behaviour. 

In particular, the second assumption implies 
that fermions both do not interact as dynami¬ 
cal degrees of freedom in the IR, and do not even 
renormalize the (purely bosonic) effective theory 
in such a way that the universality class of the 
chiral transition changes. 

Recently, the 0(4) scenario has been ques¬ 
tioned, and the chiral transition was proposed to 
be mean field behaved. 
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2. THE GROSS-NEVEU MODEL 

The reliability of the arguments against or in 
favour of the classical reduction picture is resolved 
analytically in the 3 d Gross-Neveu (GN) model at 
finite temperature jij . Here we outline the meth¬ 
ods and the results. 

The GN model we discuss is a four-point in¬ 
teracting model of N two-component massless 
fermions, with partition function and action 

Z = J DipDip exp (-S gn (ip,ip)), (1) 

Sgn(ip,i>) = J (in -3^-2- (V’G 2 ) G), 
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where / = f Q dx o f R2 d 2 x. We choose hermi- 

tian yj in the representation with 707172 = *1 2 - 
The model (|lj) serves as a good example for 
studying a “chiral” transition because it reveals 
considerable similarities to massless QCD. In 3 d 
the role of chirality is replaced by parity, corre¬ 
sponding to the transformation 

i>(x) -► ip{x) -> -ip(-x). (2) 

The model is invariant under (flj). The symmetry 
would be broken by a fermionic mass term mi>i>. 
For sufficiently large coupling A, parity symmetry 
is broken spontaneously at T = 0, and it becomes 
restored at high temperature by a second order 
phase transition. The order parameter is given 
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by 

< ip(x)ip(x) >, (3) 


which now represents the parity condensate. 

For large N it is convenient to rewrite the 
model as a Yukawa-type interacting model, using 


exp 




da exp ( — 


1 





Doing the quadratic fermionic integration, we ob¬ 
tain 


Z = c J Va exp (— N SV(cr)), 

Sy(c r) = ~tr s \ogK(z,z)j , (4) 

with K = q-d+erl, and the trace is taken over the 
two spinor indices. We impose periodic boundary 
conditions on the auxiliary field a in temperature 
direction. 

In the form d) the GN model is studied by 
means of the saddle point expansion for large N. 
Whereas the expansion in powers of A is non- 
renormalizable, the model becomes renormaliz- 
able in the l/TV-expansion. Renormalizability 
provides a large freedom in the choice of the UV 
cutoff. For simplicity, we choose a simple momen¬ 
tum cutoff, which respects parity symmetry. 


3. PHASE TRANSITION 


The minimizing configurations a(x) m i n = fi of 
the effective action (f|) are obtained as the stable 
solutions of the gap equation 

M (lb - Jl( M 2 ,r- 1 )) =0, (5) 



d 2 q 1 
(27t) 2 q 2 + x ’ 


where T = {irTv\is G Z odd}. fi(T) is equal to 
the parity condensate (|jj) to leading order in 1/N. 


3.1. N = oo 

At zero temperature, the broken phase is real¬ 
ized for sufficiently large coupling A, and we set 
the T = 0 mass scale M by 

Ji(M 2 , oo) = < • / i(°; 00 )- 

As function of T, J\ is monotonically decreasing. 
The parity condensate p.(T) approaches zero con¬ 
tinuously at T c = M/(21n2). It is straighforward 
to show that the critical exponents defined by the 
parity condensate and its correlation functions as 
T —»■ T c at to = 0 (and m —> 0 at T = T c ) become 
7 = 1, v = 1/2, r] = 0, a = 0, (3 = 1/2 (and 5 = 3). 
These are the critical exponents as predicted by 
mean field analysis. 

3.2. FINITE N 

We solve the question of the universality class 
by applying the following two methods. 

• Dimensional reduction. For large N the 
phase transition is described by a 2d scalar 
field theory. There is more than one uni¬ 
versality class with the same symmetry. By 
computing the couplings of the effective 
model, we identify that part of the phase 
diagram that corresponds to the GN model. 

• Linked cluster expansion. High tempera¬ 
ture series expansions convergent up to the 
critical points are applied to the correla¬ 
tion functions of the effective model with 
the prescribed values of the coupling con¬ 
stants. The high order behaviour allows for 
a precise measurement of the critical expo¬ 
nents. 

Dimensional reduction. This well estab¬ 
lished technique applies to local, renormalizable 
quantum field theories in D > 3 dimensions J2j. 
For sufficiently high temperature T > To, and 
for spatial momenta k sufficiently small w.r.t. T, 
the correlation functions are given by those of a 
D — 1 dimensional field theory at zero tempera¬ 
ture, which itself is local and renormalizable. 

< 0(fei);...;0(fc„) > (6) 

=< 4>$ i); • • •; 4>{k n ) >e// +0(^, ■■■, 
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The most general action is obtained by applying 
a zero-momentum expansion to the non-static ef¬ 
fective action, 

nloc (\ — n-n,m,f(n,m) 

a eff\Yst) ~ 1 * * < i >T (3-D)/2^/ t 

[-In j T>(j) ns exp (-S{T 1/2 (j) st + cj) ns ))\, 


Table 1 

The exponents 7 and v of the GN model. 


N 

7 (A) 

2 v(N) 

46.1 

1.751(2) 

2.001(2) 

24.3 

1.752(3) 

2.001(4) 

4.68 

1.749(2) 

1.999(2) 


(this provides the loophole in the no-go theorem 
of Landsman j|]), where 'T™'™'* denotes the Tay¬ 
lor expansion w.r.t. x and y of order n and m, 
respectively, constrained by 

D — 3 

f(n, m) = (D — 1) — (—-—n + m) > 0. 


For the GN model, To {N = 00) = 0. For suf¬ 
ficiently large N, dimensional reduction works 
down to the phase transition temperature T c . We 
get with the large N expansion 


= c(T) 


+T 


N 


N 3 


where we have truncated terms of order N~ 4 . All 
couplings are positive. They are known to leading 
orders of the 1/N expansion. 

Linked cluster expansion (LCE). Accord¬ 
ing to renormalizability, we are allowed to put 
the effective model on the lattice. In a form well 
known for lattice spin models, the action becomes 
in terms of the lattice fields 

So((/>o)= ( - ( 2k ) + 

x£A 2 ' 1,2 

+ [0o + X ] Xi ^o ~ 
i =2 


LCE provides convergent series expansions of 
connected correlations and susceptibilities in 
powers of k. For fixed A i we have k — T -1 / 2 * . 
Critical exponents are obtained from the high or¬ 
der behaviour of the coefficients. For instance, 
if 


X2 =< 0o(0); 0o(O) >= ^2 a L ( 2k ) l , 

L > 0 


a critical behaviour as — ( K c — re) 7 is obtained 
from 


aL 

aL -1 


2re c 


l + ^ + ^T- 1 ) 


The coefficients are computed to the 20th order. 

Some results on the exponents for various N are 
summarized in Table 1. As a result, the parity 
restoring transition of the GN model belongs to 
the universality class of the 2 d Ising model. 


4. SUMMARY 

For the Gross-Neveu model with large number 
N of flavours, dimensional reduction works down 
to the phase transition. The effective 2 d model 
has a rich phase structure, including both Gaus¬ 
sian and Ising model behaviour. However the val¬ 
ues of the coupling constants are determined by 
the reduction technique. As a result of the com¬ 
putation, the parity transition of the Gross-Neveu 
model corresponds to the Ising part of the phase 
diagram. 

For QCD we know that dimensional reduction 
works at least down to 2 T c (QH and references 
therein). Whether it works down to T c or not is 
an open question. There are additional conver¬ 
gence factors 1/T compared to the Gross-Neveu 
model because of the higher dimension, so N is 
no more required to be large. However, the renor¬ 
malized gauge coupling constant becomes large 
by approaching T c , so that perturbation theory 
normally used in computing the effective action 
becomes worse. Other computational techniques 
are required to realize the reduction below 2 T c . 
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